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SUMMARY 

This  is  a  three-part  final  report  on  the  research  supported  by  the  Air  Force  Office  of 
Scientific  Research  between  1998  and  2001,  under  Grant  F49620-99-1-0007  entitled  Sim¬ 
ulation  of  the  Transient  Aeroelastic  Response  of  a  Realistic  Aircraft  Configu¬ 
ration  During  Three-Dimensional  High  G  Maneuvers. 

1.  Motivations  and  research  plan 

Our  long-term  objective  has  been  the  development  of  a  liigh-fideUty  and  high- 
performance  simulation  capability  for  predicting  and  optimizing  the  dynamic  aeroelastic 
response  of  a  fighter  during  three-dimensional  high-G  maneuvers  in  subsonic,  transonic, 
and  supersonic  airstrearns.  Onr  focus  has  been  on  Air  Force  problems  involving  a  modern 
fighter  or  bomber,  and  relevant  to  new  approaches  for  flutter  testing,  mitigation  of  limit- 
cycle  (LCD)  and  pilot  induced  (PIO)  oscillations,  as  well  as  performance,  optimization. 
Our  starting  point  has  been  the  unique  aeroelastic  simulation  capability  developed  at  the 
University  of  Colorado  under  the  sponsorship  of  the  Air  Force  Office  of  Scientific  Research, 
and  in  partnership  with  the  Flight  Test  Center  at  the  Edwards  Air  Force  Base. 

In  order  to  achieve  the  long  term-objective  outlined  above,  we  had  defined  the  follow¬ 
ing  milestones  that  have  <;onstituted  our  short-term  objectives  during  the  three  years  of 
funding. 

•  Develop  a  high-performance  finite  element  corotational  methodology  aimed  at  captur¬ 
ing  correctly  the  nonlinear  geometric  effects  during  the  simulation  of  complex  aircraft 
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maneuvers.  Indeed,  diudng  a  realistic  three-dimensional  high-G  maneuver,  an  aircraft 
undergoes  large  displacements  and  rotations.  Therefore,  in  that  case,  simulating  the 
dynamic  behavior  of  the  flexible  airframe  cannot  be  limited  to  a  vibrational  struc¬ 
tural  analysis  around  a  fixed  equilibrium  point,  as  usually  done  for  flutter  analysis. 
Rather,  the  entire  motion  of  the  structure  must  be  tracked,  which  require.?  introduc¬ 
ing  geometric  3tiffne.s3e3  in  the  structural  model  and  performing  a  nonlinear  ana.ly.?is. 
Such  a  nonlinear  analy.9is  increases  the  c;omputatiorial  complexity  of  an  aeroelastic 
simulation,  and  therefore  calls  for  a  fast  and  scalable  nonlinear  structural  solver. 

•  Develop  a  versatile  methodology  for  designing  the  control  laws  needed  for  implement- 
iirg  a  specific  maneuver,  and  simulating  the  corresponding  actuation  of  the  control 
surfaces  of  an  aircraft.  Before  this  research,  all  maneuver.9  have  been  simulated  either 
in  a  quasi-static  fashion,  or  by  driving  some  selected  points  of  the  aircraft.  The  former 
approach  is  inaccurate  in  general,  and  fails  to  capture  transient  effect.?  in  particular. 
The  latter  approach  is  not  reliable  for  the  prediction  of  the  stresses  and  strains  that 
develop  during  a  high-G  maneuver  in  the  airframe.  Piloting  the  maneuver  through 
intelligent  inputs  to  the  control  surface.s  is  the  only  reliable  method  for  simulating  the 
aeroelastic  behavior  of  an  aircraft  during  flight.  .However,  srrch  an  approach,  which 
needed  to  be  developed  and  investigated,  requires  the  re-engineering  of  existing  flow 
solvers  and  dynamic  mesh  algorithms  to  address  the  mathematical  and  practical  con¬ 
sequences  on  the  flow  variables  and  fluid  grid  behavior  during  the  opening  and  closing 
of  the  control  surfaces. 

•  Develop  a  fast,  robust,  and  general  purpose  mesh  motion  scheme  for  flow  problems 
with  moving  and  deforming  boundarie.?.  Previous  work  in  this  area  has  been  mostly 
based  on  particulm  instances  of  the  concept  of  a  virtual  elastodynamic  fluid  grid,  and 
targeted  and  applied  to  small-amplitude  mesh  motions.  During  the  previous  bind¬ 
ing  period,  the  PI  and  his  co-workers  had  improved  earlier  the  lineal  spring  analogy 
method  by  introducing  the  concept  of  torsional  springs,  which  unlike  any  other  pro¬ 
posed  mesh  motion  algorithm,  is  mathematically  guaranteed  to  prevent  cro.S3over3 
during  an  arbitrary  inesh  motion.  However,  even  such  a  robust  algorithm  is  not  .suit¬ 
able  for  maneuvers  because  in  that  case,  the  structure  undergoes  large  displacements 
and  large  rotations.  For  such  applications,  it  is  e.ssential  to  extract  first  the  rigid  body 
component  of  the  motion  of  the  flexible  aircraft,  apply  it  to  predict  the  global  position 
of  the  fluid  dynamic  mesh,  and  then  correct  the  instantaneous  position  of  the  fluid 
grid  by  deforming  it  according  to  the  deformational  component  of  the  motion  of  the 
flexible  aircraft.  Other  mesh  motion  related  problems  that  needed  to  be  addressed 
include  tlie  handling  of  mesh  shearing  during  the  opening  and  closing  of  control  sur¬ 
faces,  and  the  automatic  treatment  of  far-field  boundary  conditions.  All  these  issues 
must  be  investigated  and  resolved  before  high-G  three-dimensional  maneuvers  can  be 
accurately  simulated. 

•  Adapt  and  tailor  various  turbulence  models  and  wall  laws  for  aeroelastic  applications. 
More  specifically,  expand  the  mathematical  formulations  of  these  models  and  laws  as 
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well  as  their  corresponding  discretization  schemes  to  the  case  of  moving  grids,  tune 
their  various  parameters  to  the  applications  of  interest,  and  design  efficient  implemen¬ 
tations  on  high-performance  parallel  computers. 

•  Develop  an  intrinsically  parallel  fluid/structure  partition  analysis  procedm-e  for  speed¬ 
ing  up  the  simulation  of  the  aeroelastic  response  of  fighters  during  complex  maneuvers, 
Indeed,  before  this  research,  our  fluid/structure  coupled  solution  algorithm  —  as  well 
as  most  if  not  all  published  algorithms  —  has  been  inherently  sequential.  While  it 
has  allowed  for  intra-pamllelism  —  that  is,  for  parallel  computations  within  the  fluid 
and  structure  analyzers  —  it  did  not  allow  for  inter- parallelism  between  the  fluid 
and  structure  computations.  The  fluid  system  had  to  be  ui^dated  before  the  struc¬ 
tural  system  could  be  advanced.  For  level  flight  applications  where  the  structure 
remains  in  the  linear  regime,  this  is  not  a  serious  limitation,  because  the  compu¬ 
tational  cost  associated  with  the  structural  analyzer  is  negligible  compared  to  that 
of  the  flow  analyzer.  However,  as  stated  earlier,  maneuvering  applications  call  for 
nonlinear  structural  analyses  that  increase  significantly  the  computational  cost  of  an 
aeroelastic  simulation.  Therefore,  onr  objective  has  been  to  develop  a  computation¬ 
ally  efficient,  and  yet  higher-order  time-accurate  staggered  procedure  for  solving  the 
coupled  fluid/structure  equations  of  motion  that  features  inter-parallelism  in  addition 
to  intra-parallelism,  and  which  tlaerefore  allows  advancing  simultaneously  in  time  the 
structure  and  flow  state  variables  without  a  significant  loss  of  accuracy. 

•  Incorporate  the  treatment  of  thermal  effects  in  our  aeroelastic  solution  methodology. 
Indeed,  the  temperature  of  the  skin  of  an  aircraft.  Tubing  grows  with  the  free-stream 
Mach  number  Moo  follows 

Tskin  =  Too{l  + 

where  Too  is  the  firee-stream  temperature,  Pr  is  Prandtl’s  number,  and  7  is  the  specific 
heat  ratio.  Hence,  at  the  higher  supersonic  Mach  numbers,  Tskin  increases,  and  the 
thermal  effects  become  important  for  the  prediction  of  the  behavior  of  both  the  fluid 
and  the  structure. 

•  Augment  our  aeroelastic  simulation  technology  with  a  high-performance  and  high- 
accuracy  computational  framework  for  performing  the  sensitivity  analysis  of  a  given 
coupled  aeroelastic  system.  Such  a  framework  is  to  feature  analytical  derivations  of 
the  gradients  of  the.  relevant  semi-discrete  operators,  rather  than  lower-order  finite 
difference  versions,  address  both  direct  and  adjoint  solution  methods,  and  include 
efficient  staggered  procediues  for  solving  all  intermediate  coupled  systems  of  ecpiatious. 
Its  main  objective  is  to  expand  the  range  of  applications  of  our  aeroelastic  simulation 
capability  to  the  aeroelastic  —  rather  than  aerodynamic  or  structural  —  optimization 
of  fighters  based  on  nonlinear  flow  and  structural  mechanics  theories. 
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2.  Achievements 

During  the  three-year  funding  period,  we  have  made  significant  accomplishments  in 
the  following  areas  pertaining  to  our  short-  and  long-term  objectives  described  above. 

2.1.  A  corotational  finite  element  method  for  aircraft  maneuvering  applications 

We  have  developed  a  systematic  approach  to  the  element-independent  corotational  dy¬ 
namics  of  finite  elements  that  eases  the  implementation  and  execution  of  the  geometrically 
nonlinear  analysis  of  aircraft  structures.  We  have  succes.sfully  applied  our  formulation  to 
the  simulation  of  various  pull  and  roll  maneuvers  of  complete  fighter  configmations  where 
the  spars,  ribs,  skin,  fuselage,  hinges,  control  surfaces,  and  discrete  masses  are  modeled  by 
.suitable  finite  elements.  This  achievement  is  documented  in  our  publications  [6,19,23,27]. 

2.2.  A  two-scale  mesh  motion  algorithm  for  simulating  aircraft  maneuvers 

Torsional  springs  provide  great  robustness  for  flutter  and/or  aeroelastic  computations 
where  the  fluid  mesh  undergoes  large  deformation.^.  However,  these  springs  and  most  if 
not  all  other  elasticity  based  mesh  motion  schemes  are  neither  sufficiently  reliable  nor  suf¬ 
ficiently  performant  when  l;he  structure  undergoes  large  displacements  and  rotations,  as  in 
maneuvering.  For  this  ree.son,  we  have  also  designed  a  new  corotational-like  mesh  motion 
strategy  where  the  motion  of  the  surface  of  the  structure  is  first  decomposed  into  a  rigid 
body  component  and  a  deformatioual  one.  First,  the  rigid  body  component  is  transferred 
to  the  fluid  me.sh  using  simple  translations  and  rotations.  Then,  the  defer mational  com¬ 
ponent  is  applied  as  a  boimdary  condition  to  the  fluid  mesh  system,  which  is  then  relaxed 
to  achieve  equilibrium  using  the  torsional  springs.  Tliis  strategy  has  enabled  the  simula.- 
hioii  of  high-angle  pitching  and  rolling  of  complete  flexible  fighter  configurations,  and  has 
speeded-up  simpler  simulations  that  were  possible  using  classical  techniques  by  a  factor 
ranging  between  2  and  10,  depending  on  some  key  configuration  parameters.  It  has  been 
summarized  in  our  publications  [6,19,23,25]. 

2.3.  Further  mathematical  development  of  the  Discrete  Geometric  Conserva¬ 
tion  Laws 

We  have  established  a  firm  theoretical  basis  for  the  enforcement  of  Discrete  Geometric 
Conservation  Laws  (D-GCL.$)  while  solving  flow  problems  with  moving  meshes.  The  GCL 
condition  governs  the  geometric  parameters  of  a  given  numerical  solution  method,  and 
requires  that  these  be  computed  so  that  the  numerical  procedure  reproduces  exactly  a 
constant  solution.  Previously,  we  have  shown  that  this  requirement  corresponds  to  a  time- 
accuracy  condition.  More  specifically,  we  have  proved  that  satisfying  an  appropriate  D- 
GCL  is  a  sufficient  condition  for  a  numerical  scheme  to  be  at  least  first-order  time-accurate 
on  moving  meshes.  While  this  results  sheds  some  light  on  the  theoretical  status  of  the  D- 
GCL,  it  does  not  fully  explain  why  it  has  been  reported  in  the  literature  that  violating 
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the  D-GCL  introduces  a  weak  instability  in  the  numerical  solution  on  moving  grids  of 
Euler  flows.  We  have  also  proved  that  satisfying  an  appropriate  D-GCL  is  a  necessary  and 
sufficient  condition  for  a  scheme  to  be  unconditionally  stable  on  moving  grids.  This  new 
result  sheds  a  new  light  on  the  relationship  between  the  D-GCL  and  numerical  scalability, 
and  provides  new  means  for  predicting  the  behavior  of  a  numei  ical  scheme  when  it  does  not 
satisfy  its  corresponding  D-GCL.  We  have  documented,  this  achievement  in  our  publications 
[2,10,14,18,24]. 

2.4.  A  third-order  energy-accurate  partitioned  procedure  for  the  fully  parallel 
solution  of  coupled  fluid-structure  interaction  problems 

We  have  developed  a  mathematical  framework  for  assessing  some  important  numer¬ 
ical  properties  of  the  partitioned  procedme  chosen  for  solving  a  coupled  fluid/structure 
sy.stem  of  equations,  and  predicting  its  performance  for  realistic  applications.  Our  analysis 
framework  is  based  on  the  estimation  of  the  energy  that  is  artificially  introduced  at  the 
fiuid/structure  interface  by  the  staggering  process  that  is  inherent  to  most  partitioned  so¬ 
lution  methods.  This  framework  also  sngge.st.s  alternative  approaches  for  time-discretizing 
the  transfer  of  aerodynamic  data  from  the  fluid  .subsystem  to  the  structure  sub.system 
that  improves  the  accuracy  and  stability  properties  of  the  underlying  partitioned  method. 
We  have  applied  this  framework  to  the  analysis  of  several  partitioned  procedures  that 
have  been  previously  proposed  for  the  solution  of  nonlinear  transient  aeroelastic  problems. 
Using  two-  and  three-dimensional,  transonic  and  supersonic,  wing  and  panel  aeroelastic 
applications,  we  have  validated  this  framework  and  highlight  its  impact  on  the  design  and 
selection  of  a  staggering  algorithm  for  the  .solution  of  coupled  fiuid/structure  equations. 
Most  importantly,  we  have  also  exploited  the  mathematical  features  of  this  framework  to 
design  a  new,  third-order  energy-accurate,  partitioned  procedure  for  the  solution  of  non¬ 
linear  transient  aeroelastic  problems  that  features  both  inter-  and  intra-parallelism.  For 
maneuvering  applications  where  the  structure  requires  a  geometrically  nonlinear  analy¬ 
sis,  this  new  algorithm  which  allows  advancing  simultaneously  in  time  the  structure  and 
flow  state  variables  has  improved  computational  speed  by  a  factor  two,  as  expected.  This 
achievement  is  documented  in  our  publications  [8,11,12,16,21]. 

2.5.  A  four-field  formulation  for  the  solution  of  aerothermoelastic  problems 

The  main  computational  difficulty  in  incorporating  thermal  effects  in  our  nonlinear 
aeroelastic  solution  methodology  has  been  the  accurate  and  efficient  solution  of  the  four 
coupled  x>artial  differential  equations  that  govern  the  flviid,  the  structure,  the  heat  transfer 
in  the  structure,  and  the  mesh  motion.  To  address  this  challenging  problem,  we  have 
extended  our  three-way  coupled  formulation  of  fluid/structure  interaction  problems  to  a 
four-field  formulation:  the  fluid,  the  structure,  the  fluid  dynamic  mesh,  and  heat  transfer. 
In  this  new  approach,  the  wall  boundary  of  the  flow  is  no  longer  assumed  to  be  isothermal 
or  adiabatic.  Rather,  the  temperature  distribution  on  the  surface  of  the  structure  is  con¬ 
sidered  as  an  unknown  that  is  determined  by  enforcing  the  continuity  of  the  temperature 
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field  and  t,eraperatm-e  fluxes  at,  the  fluid /structure  interface  boundary.  For  two-dimensional 
problems,  we  have  developed  the  semi-discrete  equations  governing  this  four-field  problem, 
and  constructed  several  partitioned  procedures  for  time-integrating  them.  We  have  applied 
this  four-field  formulation  to  the  aerothernioelastic  stability  analysis  of  flat  panels  in  su¬ 
personic  turbulent  flows,  and  to  the  thermal  analysis  of  an  F-16  multicell  wing  section 
partially  filled  with  kerozene.  We  have  documented  this  achievement  in  our  publication 
[26]. 


Time 

Influence  of  thermal  effects  on  aeroelastic  loads 
(F-16  typical  wing  section) 


T-  I  I  I  - ^ 

F-16  multicell  wing  .section 
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Time  evolution  of  the  temperature  distribution  in  the  F-16  skin 
(taking  into  account  the  presence  of  fuel) 


2.6.  Analytically  based  sensitivity  analysis  and  optimization  of  nonlinear  aeroe- 
lastic  systems 

We  have  considered  the  problem  of  optimizing  for  steady-state  conditions  a  given 
aeroelastic  system  by  varying  both  aerodynamic  and  structural  parameters  such  as  the 
shape  of  the  dry  or  wet  surface,  and  the  orientation  of  the  composite  fibers.  We  have 
developed  a  solution  methodology  for  this  optimization  problem  that  is  based  on  a  sequen¬ 
tial  quadratic:  programming  procedure  where  the  gradients  of  the  optimization  criteria 
with  respect  to  the  optimization  variables  are  determined  by  an  analytical  approach.  In 
this  methodology  which  features  both  the  direct  and  adjoint  approaches  for  optimization, 
we  evaluate  at  each  optimization  step  the  aeroelastic  steady-state  response  of  the  system 
using  a  staggered  procedure  that  couples  efficiently  a  finite  element  solution  method  for 
the  structure  subsystem  and  a  3D  Euler  finite  volume  method  for  the  fluid  subsystem.  We 
have  illustrated  this  optimization  methodology  with  several  three-dimensional  examples 
and  documented  it  in  our  publications  [4,5,17,20]. 

2.7.  Simulation  of  the  aeroelastic  behavior  of  complete  F-16  configurations  and 
validation 

We  have  apidied  our  .simulation  technology  to  the  flntter  clearance  of  an  F-16  Block 
40  in  clean  wing  configuration  but  with  tip  missiles,  for  0.7  <  Moo  <  1-4  at  the  altitude 
of  3,000  m.  Based  on  modeling  information  provided  by  Lockheed-Martin,  we  have  con¬ 
structed  at  the  University  of  Colorado  a  detailed  but  undamped  three-dimensional  FE 
structural  dynamics  model  for  the  F-16  Block  40  in  clean  wing  configuration  but  with  a 
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missile  and  launcliing  system  at  each  wing  tip.  This  FE  model  features  bar,  beam,  solid, 
plate,  shell,  metallic  as  well  as  composite  elements,  and  a  total  of  168,799  dels.  It  re¬ 
produces  correctly  the  first  ground  bending  and  tor.sion  frequencies  which  were  measured 
as  4.76  Hz  and  7.43  Hz,  respectively.  Using  F-16  CAD  data  provided  by  the  Air  Force 
R.esearch  Laboratory  at  Wright  Patterson  and  ignoring  the  wing  tip  missiles,  we  have  gen¬ 
erated  a  surface  grid  with  63,044  grid  points  then  an  Euler  fluid  volume  rne.sh  with  403,919 
vertices.  To  determine  the  bending  and  torsional  aeroelastic  parameters,  we  have  em¬ 
ployed  a  procedure  replicating  what  occurs  in  flight  testing.  We  have  excited  the  structure 
in  an  appropriate  manner  and  simulated  its  re.sponse  to  the  prescribed  initial  disturbance. 
For  each  different  Mach  number,  thi.s  generated  168,799  signals,  one  for  each  dof  of  the 
detailed  FE  structural  model.  We  have  computed  these  signals  for  a  little  more  than  two 
(;ycle.s,  then  applied  the  Eigensystem  Realization  Algorithm  to  axtract  from  a  few  of  them 
representing  sensor  information  •  the  sought-after  aeroelastic  pa.ra.meters.  We  have 
compared  our  result.s  to  flight  te,st  data  provided  by  the  Edward  Air  Force  Base  and  ob¬ 
served  less  than  7%  relative  errors.  This  very  good  correlation  has  attracted  the  attention 
of  major  aerospace  companies  with  whom  we  are  now  working  on  further  validations. 


Detailed  finite  element  structural  model  of  an  F-16  Block  40 
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